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ABSTRACT 



’"our satellite esti nation techniques were applied to 
tv;elve cases durinp the v/inter/sprinj (19'")6-S7) alonq the 
v/Gst coast of the "’nited States. These estinates v/ere 
then compared with observed precipitation patterns fron 
radar and s’lrface observations to evaluate their skill. 

esults proved the National En vironnental Satellite, 
kata and Information Service (NESDIS) subjective 
technique superior. Other methods v;ere Ik threshold, 
color bispectral and the h'aval Postgraduate School (?IPS) 
bispectral minicomputer program. All techniques 
overestimated observed precipitation, but the NESDIS 
method substantially reduced the amount of over- 
estimation. All techniques showed excellent alignment 
and \s'ere more successful representing the western edge 
rather than the eastern edge of precipitation. Only the 
;!ESDIS and 'IPS bispectral approach showed skill in 
depicting precipitation intensity. 



A 



TABLE OF COL’TENTS 



I. INTRODUCTION 11 

II. EONITORING RAINFALL: THE NATURE OF THE 

y’RO^LE:: 15 

A. RAINFALL NEASUREMENT 15 

B. NEED OF SATELLITE :!0::iT0RING 17 

III. ENPERINENTAL r’ROCEDURE 19 

A. EVALUATION PROCEDURE 19 

D. TYPES OF VERIFICATION 20 

1. Radar Sumnary Chart. 20 

2. Surface Observations 21 

r: . EXPERIMENTAL TECHNIQUES 22 

A. NESDIS TECHN'IOUE FOR SHORT RANGE 

FORECASTING OF PRECIPITATION 2 2 

B. SATELLITE DATA PROCESSING AND DISPLAY 

SYSTEM (SPADS) TECHNIQUES 23 

1. Infrared Thresholding 23 

2. NPS Geostationary Satellite Cloud 

and Precipitation Analysis Model. ..25 

3 . Bispectral 25 

V. EXPERIMENT RESULTS 27 

A. CASE NUMBER 9 (EXCELLENT AGREEMENT) 30 

B. CASE NUMBER 1 (POOR AGREEMENT) 3 5 

C. SUMMARY OF ALL CASES 39 

•'I. CONCLUSIONS AND RECOMMENDATIONS 43 

APPENDIX A; TABLES 45 



5 



APPENDIX B: FIGURES 



59 



LIST OF REFERENCES 

INITIAL DISTRIBUTION LIST 79 



6 



LIST OF TABLES 



1. Comparison of Techniques for Case Number 1 

at 1830 G':T on 4 Pecember 1986 43 

2. Comparison of Techniques for Case Number 2 

at 1830 GMT on 22 December 1986 44 

3. Comparison of Techniques for Case Number 3 

at 1830 GMT on 23 December 1986 45 

4. Comparison of Techniques for Case Number 4 

at 1830 GMT on 1 January 1987 46 

5. Comparison of Techniques for Case Number 5 

at 1830 GMT on 2 January 1987 47 

6. Comparison of Techniques for Case Number 6 

at 1830 GMT on 3 January 1987 48 

7. Comparison of Techniques for Case Number 7 

at 1800 GMT on 22 January 1987 49 

8. Comparison of Techniques for Case Number 8 

at 1830 GMT on 2 February 1987 50 

9. Comparison of Techniques for Case Number 9 

at 2130 GMT on 9 February 1987 51 

10. Comparison of Techniques for Case Number 10 

at 2130 GMT on 12 February 1987 52 

11. Comparison of Techniques for Case Number 11 

at 1830 GMT on 13 February 1987 53 

12. Comparison of Techniques for Case Number 12 

at 0030 GMT on 13 March 1987 54 

13. Summary of Tables 1-12 55 

14. NESDIS Category and Subjective Skill Compared 

with Radar and Surface Observations 56 



7 



LIST OF FIGURES 



1. GOES Visible ;iercator Image for 2130 G'AT on 

February 1987 57 

2. GOES Infrared Mercator Image for 2130 G!!T on 

9 February 193 7 58 

3. EMC Radar Summary Chart for 2135 GMT on 

9 February 1937 59 

4. Infrared Schematic of the NESDIS Cloud Band 

Category 60 

5. Infrared Schematic of Evolution - Subsynoptic- 
Scale L’ave category - ^vaves in an Overrunning 

Zone 61 

6. NESDIS Technique Overlayed on High Resolution 

GOES-6 Visible Image for 2100 GMT on 9 February 
1987 62 

7. Infrared Threshold (-25^ C to -55^ C) Image 

for 2130 GMT on 9 February 1987 63 

8. NFS Model Results for 2130 GMT on 9 February 

1987 64 

9. Bispectral Image for 2130 GMT on 9 Februarv 

1987 1 65 

10. GOES Visible Mercator Image for 1830 GMT on 

4 December 1986 66 

11. GOES Infrared Mercator Image for 1830 GMT on 

4 December 1986 67 

12. NMC Radar Summary Chart for 1830 G!1T on 4 

December 1986 68 

13. NESDIS Technique overlayed on High Resolution 

GOES-6 Visible Image for 1901 GMT on 4 December 
1986 69 



8 



14. Infrared Threshold (-25^ C to -55*^ C) for 

1830 GMT on 4 Dece-nber 1986 70 

15. '-'PS Model Results for 1330 GMT on 4 December 

1987 71 

16. Bispectral Image for 1830 GMT on 4 December 

193 7 7 2 

17. Comparison of Area Coverage by Technique 73 

13. Cloud-Top Temperatures Correlating with 

"Observed Precipitation 74 



9 



A C !< ' ' 0 1 •/ L K P G F ' 1 F ' I T 



I would like to express ny appreciation to ny wife, 
Faye anci :~:y daughters, Shelley and ."^age, for supporting 
'le through the long nights away from them spent 
collecting data and \/riting. I would also like to thank 
Dr. Carlyle h'ash, ny thesis advisor, for his insight, 
knowledge, and guidance and Dr. Philip Durkee for acting 
as my second reader. In addition, I \7ould like to express 
ny gratitude to /'r. Lang Chou for spending many hours 
reprogramming the 1^’aval Postgraduate School Geostationary 
Satellite Cloud and Precipitation Analysis '!odel so it 
would run on the computers at the !Iaval Environmental 
Prediction Research Facility (NEPRF). And finally, I 
would like to acknowledge the :1EPRF Command and 
especially Dr. Andy Goroch, Mr. Ron Picard and Mr. Jim 
Clark for their support and assistance. 



10 



I'lT^.ODUCTIOri 



I . 

For the past decade, Naval Oceanography was 
centralized around Fleet '^umerical Oceanography Center 
(FNOC), llonterey, CA. Alnost every analysis or forecast 
oroduct was derived or originated from FNOC* The volume 
of products became so large that their computer facility 
had little capacity for future requirements. Thus, the 
current goal of the xNaval Oceanography community has been 
to diversify computer assets by allowing the oceanography 
centers to maintain computers capable of modifying FNOC 
products and/or generating their own regional products. 

In addition to providing the user with more guidance, it 
also aids the survivability of the entire oceanography 
netv;ork in case of wartime conditions. 

As a result of this general diversification, 'Javal 
Oceanography regional processing capabilities will 
rapidly improve over the next five years. This will be 
accomplished through the development, procurement and 
distribution of three major equipment suites to the 
fleet. They are: the AN/SMO-11 Satellite Receiver, the 
N'aval Environmental Data Network (NEDM) Oceanographic 
Data Distribution and Expansion System (NODDES) / 



Satellite-Data Processing and Display System (SPADS), and 
the Tactical Environmental Support System (TESS). All 
three systems have in common satellite processing and 
display, \vhich enables the oceanographer on the ship or 
at the remote station to make evaluations of current 
weather parameters wit ii out the assistance of regional 
centers (e.g. Guam, Pearl Harbor, Norfolk or Rota) or 
F'lOC. In addition, all three systems are scheduled to 
arrive in the fleet almost simultaneously. 

Although ship, aircraft and weapon system technology 
has advanced steadily over the past decade, oceanography 
and meteorology continue to play a primary role in the 
planning and execution of major operations. Satellite 
capabilities provide the operational command with a 
passive (i.e., the ability to maintain communication 
silence) means of obtaining numerous oceanographic and 
meteorological parameters (e.g., cloud, temperature, 

\v'ind, moisture, pressure, sea condition, sea-surface 
temperature, weather element, etc.). The accuracy and 
extent of these observations have improved dramatically 
with better satellites, sensors and applications. Thus, 
it becomes imperative that the Naval Oceanography 
community utilizes the latest techniques for application 
of satellite data. 
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This thesis is designed primarily to evaluate 
satellite processing and display capabilities for the 
detection, identification and classification of 
precipitation over the west coast of the United States. 
The accuracy of four satellite precipitation estimation 
approaches will be evaluated using twelve win ter /spring 
cases. Three methods use the Satellite-Data Processing 
and Display System (SPADS) software pad: ages v;hile a 
fourth is primarily subjective. Verification is conducted 
with the National Meteorological Center^s (^N^C’s) radar 
summary charts and surface observations. 

Chapter II of this thesis reviews the problem of 
monitoring rainfall including key research involving 
satellite derived precipitation. Chapter III gives a 
detailed description of the experimental procedure. 
Chapter IV includes the four precipitation estimation 
schemes including the National Environmental Satellite, 
Data and Information Service (NESDIS) technique for 
identifying precipitation, SPADS bispectral and 
thresholding techniques, and the current Naval 
Postgraduate School (NPS) cloud and precipitation model. 
Chapter V describes the experimental results using a 
case demonstrating excellent agreement and a case 
demonstrating poor agreement. The experimental results 
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for these two cases are provided in detail. The remaining 
ten cases, analyzed with the sane attention to detail, 
have been summarized in tables for ease in reading. 
Conclusions and recommendations for all twelve cases are 
found in Chapter VI. 
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II, 



MONITORING RAINFALL: 



t^it? 



OF th: 



A, RAINFALL MEASUREMFNT 



"^arrett and Martin (1981) describe precipitation as 
one of the most variable elements in v/eather. It varies 
with respect to its frequency, duration, intensity and 
spatial pattern, not to mention its appearance in both 
solid and liquid forms. Rainfall organization and 
distribution, influenced by atmospheric circulation, are 
greatly enriched by global geography. These factors 
result in rainfall distribution patterns of great 
complexity. 

Conventional techniques for measuring precipitation 
include the traditional ^ situ devices (rainguages) and 
ground-based remote sensing systems (weather radar), both 
of which have been used successfully during the past 
twenty-five years. Although these methods provide 
reliable rainfall data, it is not surprising that 
monitoring with the desired detail, timeliness and 
accuracy is often difficult. 

Measurement of rainfall by gauges depends on 
topography, site, wind and gauge design. In addition, 
accessibility of locations, availability of personnel to 
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read the instrument, suitable pov/er supply, security of 
station, ability to communicate readings quickly and 
accurately, and lack of observation platforms over the 
oceans magnify the difficulty of using guages to 
observe precipitation. It is easy to see why better 
methods are sought. 

In the case of the remote weather radar, problems 
with resolving the proper relationship of back scattered 
microwave energy to drop spectrum size, attenuation, 
absorption, reflection and signal calibration all 
present difficulties in the measurement of precipitation 
intensity and amounts. Barrett and Martin (1981) indicate 
that spatial correlation of rainfall decreases rapidly 
with increasing distance between points (e.g., rainfall 
at points separated by distances as small as 50 km is 
practically unrelated). Further discussions of the 
capabilties and weaknesses of radar measurement are found 
in Mason (1971), Battan (1973), and Hudlow et al. (1979). 
Hov/ever, Barrett and Martin point out that radar does 
b.ave an advantage over conventional gauges. It is 
spatially continuous over its prescribed range, usually 
125 miles over flat terrain. 
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NEED OF SATELLITE MONITORING 

Precipitation plays an inportant role in weather, 
climate, hydrology and oceanography. As a primary part 
of the hydrological cycle, man depends on it as a natural 
source of generally pure \/ater for use in agriculture, 
forestry and commerce, v;hich includes transportation, 
communications and industry. It also supplies the net of 
heating from condensation within the atmosphere, \'/hich is 
an important input to any global atmospheric/oceanic 
model. Thus, the significance of monitoring this weather 
parameter globally grows as modeling techniques improve. 

The Navy’s interest in precipitation mirrors many of 
the reasons mentioned above, but also includes tactical 
requirements. Unfortunately, conventional methods of 
collecting rainfall data provide little or no coverage 
over the oceans. Land coverage, often sparse at best, may 
disappear during a confrontation, such as the recent 
Falkland Island conflict v/here Argentina stopped 
reporting surface observations for the duration of th. e 
conflict with England. Other tactical decisions (i.e., 
ship movements, amphibious landings, carrier operations, 
etc.) require current information on cloud cover, 
precipitation, winds and temperature as v/ell as the 
ability to make accurate short and long range forecasts. 
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Therefore, the !'avy must strive to upgrade its ability to 
monitor these parameters, especially over the oceans. 

Other reasons for using satellites vice conventional 
methods include the ability to resolve larger spatial 
patterns, timeliness of observations and a global 
data base, h'ith the new '"aval Oceanography systems, 
digital satellite data v;ill soon be available to the 
remote user. Satellite methods for estimation of 
precipitation will enhance the final analysis and 
forecast products for both regional and local sites. 

This thesis will focus on comparing satellite 
precipitation schemes which provide reliable analysis of 
synoptic precipitation patterns. 
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Ill . 



^AL P’lOCEDU^ 



A . PVALUATIO?: PROCEDURE 

The central v/est coast of the Pnited States served 
as the evaluation region for this investigation. This 
geographical location v/as selected for two reasons. 

Pirst, winter maritime weather systems, only slightly 
modified by topography, move through this region on a 
rather routine basis. And second, two National Weather 
Service (NWS) radar observation stations are located 
near the center of this area to provide reliable 
verification of satellite imagery through use of 
automated radar summary charts. 

The satellite methods of mapping precipitation 
studied were based on visible and infrared imagery. 
Although passive microwave techniques are being 
developed, there were no operational microwave data 
available. Futhermore, the GOES imagery could be 
collected every 30 minutes v;hich meant images could be 
selected which closely coincided with Radar Summary Chart 
times. This v;ould not have been possible with microwave 
data since they come from a polar orbiting satellite and 
orbit times vary daily. 
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GOES visible and infrared imagery was collected on 
the Global Acquisition and Data Handling System (GADIIS) 
at the Naval Environmental Prediction Research Facility 
(NEPRF) in Honterey, CA* The data v;ere collected over a 
five-month period from November 19S6 to !'.arch 1937. From 
that data set, tv;elve cases were selected by the 
following decision criteria. If automated radar summary 
charts indicated precipitation \;ithin the region under 
investigation, GOES imagery along v;ith temperature field 
data from FNOC were collected. If the imagery and field 
data were collected with no problems, then surface data 
\^ere ordered for the closest synoptic hour. A satellite 
analysis of the imagery was completed using four 
techniques. These techniques are described in the next 
chapter. Finally, comparisons between results using 
these techniques and actual radar and surface 
observations were made and evaluation statistics derived. 

3. TYPES OF VERIFICATION 

1 . Radar Summary Charts 

The automated radar summary chart, produced by 
the National Weather Service, was a major verification 
tool for this project. Imagery were received to coincide 
as near as possible with the times on the radar summary 
chart. The radar information was then transferred 
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manually from the radar summary chart, on a polar stereo- 
graphic projection, onto a GADIIS mercator background 
of the eastern north Pacific ocean and the v/estern United 
States. A viewgraph of the precipitation pattern \/as 
talcen and compared v/ith satellite estimation techniques. 

2 . Surface Observations 

Surface observations were used to supplement t!ie 
radar summary charts. Data for the nearest synoptic hour 
v/ere collected from UPS’s archives and FNOC’s data hanlc. 
UPS archives provided surface charts with excellent 
coverage over the area of interest. FNOC charts provided 
sparse data over the area of interest but added ship 
observations to the immediate west. Together, the two 
sources provided good coverage of the area of interest 
and adjacent ocean. For some cases, however, surface 
observations were given more weight than the radar 
summary chart since the images and surface observations 
occurred at the same exact time. In these cases the radar 
summary chart was valid 30 minutes following the image 
time. 
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IV. 



TA?E^.i::T:.TkL tecii::ioues 
A. ::esdis technique for short ra::ge forecasting of 

?RECI?ITATIO:i 

The NE3DIS technique, developed by Scofield and Spayd 
(1934), uses geostationary satellite inagery for 
estimating precipitation from extratropical cyclones 
during the v; inter season. They divide extratropical 
cyclones into five categories. Each category is defined 
by distinct cloud pattern evolution and IR enhancements 
■/hich indicate where the light, moderate and heavy 
precipitation should occur. From this analysis, short 
range (1- to 3-hour) forecasts are prepared. These 
techniques are very subjective and usually overestimate 
precipitation intensities and amounts, as indicated in 
Scofield et al. (1932) and Scofield and Spayd (1983). 

Scofield and Spayd (1934) schematics of cloud 
patterns associated with precipitation were developed by 
examining IR data for various types of v; inter storms. The 
schematics are used as aids in analyzing the location and 
magnitude of continuous or sho\:ery precipitation using IR 
threshold and cloud pattern structure \vrithin the v/ inter 
storm clouds. Images of this technique, produced 
operationally by the regional National Weather Service 
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(NUC) forecaster at Redwood City, CA, and sent out over 
the UNIFAX circuit, will be compared to observed 
precipitation patterns in the tv/elve validation cases. 
There are five basic types of cloud patterns. "I’hey 
are: connca head category, baroclinic leaf category, 
subsynoptic-scale wave category, cloud band category, and 
overruning category. More details on these categories are 
provided in the case by case analysis. The MESDIS method 
differs substantially from the objective methods in that 
it is a product with human input. The comparison of both 
subjective and objective approaches v/ill assess the role 
of human decisions in satellite precipitation estimation, 

3. SATELLITE DATA PROCESSING AND DISPLAY SYSTEM (SPADS) 
TECIi::iQUES 

1 . Infrared Thresholding 

Infrared thresholding is the original subjective 
technique (i.e., digital data transformed to shades of 
gray or colors) used for mapping precipitation, however, 
due to its inability to easily distinguish between thin 
and thick cold clouds, Barrett and Martin (1981) 
conclude that IR thresholding is only effective v/hen 
v/orking with isolated convective activity. Otherv/ise the 
technique consistently overestimates the observed 
precipitation pattern. 
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There are nunerous inaf^e enhancements designed to 
isolate precipitating clouds. Current objective methods 
use r^ore sophisticated means of thresholding to 
’istinguisli types and intensity of precipitation, because 
all of the methods described here use IR data in 
threshold form, it is presented here as an analysis 
technique . 

A range of IR thresholds are used in identifying 
a rain/no rain condition. Some models, like the NFS 
model, use two threshold cutoffs (i.e., 247.5 K and 223 
K) for discriminating light and moderate intensities of 
precipitation, ^ue to the many accepted standard 
thresholds, this technique produces an image of a 
specific range of cloud-top temperatures (i.e., 218 K to 

25 <S K) likely associated with precipitation. The range 
v;as selected after reviewing articles by Scofield and 
Oliver (1977), Megri and Adler (1981), Paul (1983), and 
Spray (1935) which documented their IR threshold ranges 
of 193 K to 243 K, 229 K to 260 K, less than 247.5 K 
to 254 K and less than 246 K to 254 K respectively. 

In addition, a range of cloud-top temperatures 
are recorded for each case study to be presenterd in 
figure form later in the thesis. The ranges coincide 
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v/ith active precipitation areas fron the radar summary 
chart and surface observations. Then, these ranges are 
compiled and presented in graphical format for use in 
comparing past east coast results with this season^s 
^/est coast results. 

2 . ''PS Geostationary Satellite Cloud and 
r r e c i p 1 1 a t i o n Analysis i ! o ! e 1 

T’ash et el. (198 5) describes an automated cloud 
and precipitation analysis scheme. The precipitation 
estimation method is a modified version of Liljas*s 
(1932) threshold technique adopted from the results of 
Muench and Teegan (1979). The precipitation module 
activates if the cloud-type registers as nimbostratus or 
cumulonimbus. Intensity depends on infrared and visual 
values; the colder infrared temperatures and brighter 
visual albedos produce heavier intensities. As indicated 
before, the infrared precipitation threshold consists of 
two values (247.5 T and 223 K), v/hich equate to light and 
moderate precipitation respectively. The values vary 
slightly because they are dependent upon observed 
temperature profiles. 

3 . Bispectral 

Due to time constraints, it may not be possible 
at sea using a TESS system to execute a complete cloud 
analysis model. A simpler bispectral approach, using both 
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visible and infrared data, has been developed for SPADS 
software* The bispectral technique combines both visible 
and infrared images and overlays the two images in 
different colors (e.g., IR-red and VTS-creen). ^'here t!ie 
images cor-ibine to produce different shades of orange and 
y ell o\>: , the technique indicates significant cloud 
thiclcness (visual) and cold (high) cloud-tops (infrared), 
givin^ a strong probability of precipitation. The closer 
the shaciing comes to a pure yellow the greater the 
thickness and the colder the cloud-top temperature, which 
’would indicate a heavier intensity. 

The SPADS high quality color display is an 
efficient method of showing bispectral satellite data. 
Darrett and !^[artin (1931) discuss key precipitation 
studies using bispectral methods. The most successful 
application has been by Lovejoy and Austin ( 1979), v/here 
radar d. ata guides the choice of critical thresholds. 
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V. 



"^XPr^.I’:E:^TAL RESULTS 



The evaluation results are now presented. To 
illustrate the evaluation procedure, attention is 
'focused on tv;o of the t\/elve cases, a case with 
excellent and a case v;ith poor agreenent. Results from 
all twelve cases are then sunr.arized in Tables 1 - 12. 
Each case looks at all four of the experimental 
techniques mentioned in the previous chapter and makes a 
subjective comparison of actual radar and surface 
observation coverage. 

The satellite estimates of precipitation from each of 
the four techniques is compared to the observed radar 
data for the central west coast region of the United 
Rtates. An overlay of the radar coverage is placed over 
the satellite estimates. Attention is first focused on 
the area of the satellite estimate that is verified by 
the radar and surface data. Preferably, satellite 
estimation techniques should have been compared with 
digital radar observations to generate objective threat 
scores (i.e., comparison of two fields or patterns; 
amount of overlapping area divided by the overall area of 
both fields to yeid a percentage). !Iowever, digital radar 
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observations could not be obtained to correspond to the 
inage times. Therefore, a subjective threat score is 
determined by comparing the correct area v/ith the overall 
analysis and verification areas. T\/o types of patterns 
\;Qre found in these evaluations, '^he majority of 
satellite estimates were larger than the observed radar 
or surface observation pattern. So the threat score was a 
percentage of the satellite estimate. In the few cases 
where the radar and surface observation areas were larger 
than the satellite estimates, the threat score was the 
percent of the rain area correctly classified by the 
satellite method. The percentage is evaluated 
subjectively v/ith an estimation error of approximately 
ten percent. Since the percentage of agreement is 
subjective, four categories are utilized to describe the 
results. The categories are: Excellent Agreement -greater 
than 75%, Good Agreement -betv/een 50-75%, Fair Agreement 
- between 25-50%, and Poor Agreement - less than 25%. 
Text, each technique is compared with the observed 
precipitation and classified as either an over- or 
underestimation. In tv/o cases (i.e., 2 and 11), it 
overestimates the northern portion of the frontal band 
u'hile it underestimates the southern tail of the frontal 
band. 
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A third evaluation criterion is the alignment of 
satellite-estimated precipitation compared to that of the 
observed area. The criterion for deciding alignment 
involves the orientation of the overall satellite 
estimates compared with the orientation of the radar and 
surface observation data. If the alignment axis of the 
precipitation pattern varied by more than 20 degrees 
(i.e., north to north-northeast), it was considered 
incorrect. Otherwise the alignment was considered 
correct . 

Precipitation intensity v/as also verified by 
subjectively comparing the satellite estimation with the 
information provided on the radar summary chart and 
surface observations. The evaluation v/as characterized as 
either an overestimate, underestimate or correct. 

A fifth evaluation criteria was the accuracy of the 
leading edge of precipitation. The leading edge of 
precipitation, normally the eastern edge, coincided with 
the onset of precipitation and is a critical forecast 
variable. The leading edge was considered correct if the 
satellite estimate and observed data agreed within two 
degrees of longitude. Otherv/ise, the technique's leading 
edge was labelled either substantially east or v;est of 
the observed precipitation edge. 
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Since radar sunmary charts v/ere initially selected 
as the nost representative method of presenting 
precipitation patterns for comparison, most of the 
imagery v;ere taken at 30 minutes after an intermediate 
synoptic or synoptic hour !i i c h \ ; a s v; i t h i n five minutes 
of the radar summary chart tine* Thus, surface 
observations, taken approximately ten minutes before the 
!iour, provide only a secondary tool in making these 
qualitative comparison* Hov/ever, there are two in which 
the imagery were taken exactly on the synoptic hour. In 
these cases more credence is given to the surface 
observation charts* 

A* CASE :!U:iBER 9 (EXCELLENT AGREEMENT) 

The Rational !Ie teorological Center's (H!IC's) Mean 
Sea Level (!I3L) pressure analysis, for 2300* GMT on 9 
Eebruary 19S7 (not shown), shows a fast moving cold front 
located along the west coast of the United States* It is 
approaching a weak wave originating from an upper-level 
perturbation located over northern California* Weak high 
pressure dominates the remainder of the western United 
States* The visible and infrared images. Figs* 1 and 2 
respectively, depict the complex cloud structure 

”!!C failed to produce a 2100Z surface analysis. 
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associated v;ith this situation. The visible inage shows 
a distinct comma cloud v/ith its center located at 
approximately 49.0^ h and 131.5^ W. Toth visible and 
infrared images indicate an unusual bulge in the frontal 
cloudiness over nortliern California. 

The radar summary chart for 2135Z (Fig. 3) shov/s 
light to moderate rain and rainshowers occurring from 
Astoria, 03 southward to Los Angeles, CA. The entire area 
of precipitation is m.oving northeast at 25 to 30 knots. 
'Rainfall intensities increase over northern California, 
probably in relation to the v;eak v/ave developing in that 
region. Coastal precipitation observations, depicting 
rain and rainshowers from Astoria, OR southv/ard to 
Los Angeles, CA, increase confidence in the area 
precipitation coverage portrayed by the radar summary 
chart . 

The CTSDIS technique, developed by Scofield and Spayd 
(1984), identifies the frontal zone as a Cloud Band 
pattern. The IR schematic of the RESDIS Cloud Rand 
Category can be seen in Fig. 4. The weak wave was 
identified as a Subsynoptic-Scale ".'ave in an overrunning 
zone. Its IR schematic can be seen in Fig. 5. The 
schematics of evolution of cloud patterns associated v;ith 
moderate to heavy precipitation are used by the IL.7S 
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Redwood City forecaster to aid in analyzing the location 
and magnitude of the continuous or showery precipitation, 
Scofield and Spa yd (1984) describe both categories as 
f ol lov;s : 

Cloud Rand: Cirrus clouds comprise most of the cloud 
band, h o v; e v e r , t Ii o s e bands which possess low-level 
tails are active (contain precipitation). In addition, 
\;hen newer cloud bands develop in the rear portions of 
the baroclinic zones, the older bands in the front 
portions of the baroclinic zones often dissipate. 

Subsynoptic-Scale b’ave: At times the v/ave appears to 
develop in an area of convection along the southern 
portion of the baroclinic zone and propagates 
nor theas tv/ard along the zone. 

Comparisons of the operational Redwood City hSSDIS 

product, as seen in Fig. 6, v/ith observed radar and 

surface observation precipitation patterns indicate 

estimates that are classified as excellent agreement. By 

the size and intensity of shower activity, the 

Subsynoptic-Scale I/ave schematic (Fig. 5) using T+4 hours 

shows the best correlation. The alignment of the 

precipitation patterns is extremely good, intensities 

showing light rain and rainshowers are slightly 

overestimated, and overlapping area coverage as related 

to threat score is close to 90%. The leading edge of the 

precipitation pattern is excellent except over a small 

section of northern California where the technique 
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underestimates the actual precipitation pattern slightly 
to the v;est. 

The IR threshold technique observed in Fig, 7 
is classified as excellent agreement. Overlapping area 
coverage is "^50, alignment and leading edge are exact, 
''ovever, intensities (i.e., colder cloud-top 
temperatures) are again overestimated along the Oregon 
coast. IR enhancements indicate that temperatures in the 
range of 225 T to 243 T best represent the pattern 
outlined in the radar summary chart for case number 9. 

As indicated :)efore, results from all cases v/ill be 
placed in graphical format for comparison with mean 
thresholds derived from previous studies. 

The hPS model results observed in Fig. 3 is 
classified as excellent agreement (i.e,, area coverage is 
75T), Again, the technique overestimates the leading edge 
of precipitation in the mountains northeast of Sacramento 
but other v/ise the rest of the pattern is correct. And as 
stated above, radar and surface observation coverage in 
this region is suspect. The model indicates areas \;here 
more intense precipitation is likely, and like the FESDIS 
technique, it indicates rainshowers occurring 
throughout the precipitation band. The alignment is 
correct . 
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The bispectral image observed in Fig*9 is classified 
as excellent agreement. Area coverage is 90%. Bright 
yellows highlight all areas of precipitation. However, 
the area of bright yellov/s actually exceeds the area of 
precipitation over the mountains northeast of Sacramento, 
California. Since the radar cannot see over the mountains 
and surface reports are extremely sparse in this region, 
it remains inconclusive whether the bispectral 
representation is overestimating the actual precipitation 
at all. Otherv/ise the leading edge is correct. The 
alignment of the area is exact. Precipitation intensity 
is difficult to estimate using the bispectral method 
because it is hard to determine when the shade of yellow 
is getting brighter. But it appears that the brightest 
yellows are again north of the weak wave. 

Table 9 contains a summary of all the information 
above. All methods received an excellent score, but the 
:;BSDIS and bispectral techniques scored slightly higher. 
The IIESDIS and IR techniques underestimated the area 
coverage while the other two methods overestimated it. 

All techniques overestimated the intensity of 
precipitation, shov/ed correct alignment and properly 
defined the leading edge of precipitation. 
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B. CASE ;rj;!BER 1 (POOR agree;;ent) 

The MMC MSL pressure analysis for 1800 GMT on 4 
December 19S6 (not shown) displays a weak cold front 
approaching the v;est coast of the United States. The 
visible and infrared imagery, Figs. 10 and 11, shov; a 
comma shaped cloud v/ith its center located at 40.5^ U and 
130.1^ I/. A massive cirrus shield overshoots the active 
frontal zone by nearly 700 nautical miles. A weal: 
stationary front still lingers along the border of 
k’ashington and Oregon, further enhancing the prefrontal 
cloudiness . 

The BMC radar summary chart, seen in Fig. 12, shows 
broken areas of precipitation from northern Uashington 
southv/ard to central California. The areas of 
precipitation are moving eastward at 15 to 20 knots. The 
surface analysis for 1800 GMT only indicates rainfall 
occurring at one coastal station in Oregon. Thus, it does 
not add much confidence to the radar depiction. II ov; ever, 
the radar summary chart time, 2135 GMT, is closer to the 
image time than the surface observation chart. So it is 
used as the primary verification tool in this case. 

The NESDIS technique classifies the synoptic 
situation above as a Cloud Band category. Discussed in 
case number 9 and illustrated in Fig. 4, this case 
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emphasizes the importance of cloud structure, especially 
as it relates to those bands v;ith low level tails. The 
?!ESDTS analysis technique, shown in Fig, 13, severely 
underestimates the amount of precipitation actually 
occurring along the \;est coast and therefore exhibits 
poor agreement v;ith the radar summary chart. Hov/ever, it 
does capture the second system off shore. The NESDIS 
approach in a cloud band case focuses attention at the 
immediate cloud bands (using the cloud tails as a guide). 
The large mid and high cloud shield ahead of the bands is 
correctly classified as no precipitation. The error in 
this case is to underestimate precipitation in the first 
band. Fig, 13 is a working copy of the NESDIS technique 
supplied by the NTS, Redwood City, CA forecasting center. 
The red m.arkings on the picture should be ignored since 
they represent a short range forecast and are not related 
to this study. Area coverage or threat score in the 
region of interest is approximately 15% v/hich classifies 
the estimation technique as poor agreement. The alignment 
of the precipitation pattern is also poor. However, the 
misalignment, in part, is probably due to the initial 
analysis. Although the first baroclinic zone is severely 
underestimated, there is enough depicted to clearly 
identify the leading edge of precipitation. Since a 
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stron;^,er baroclinic zone is developing to the rear of the 
first zone, intensities are correctly depicted by sho\/ing 
light precipitation occurring along the eastern 
baroclinic zone. 

In 14, the I?, threshold technique is classified 

as poor agreenent. Overlapping area coverage is 20h'. The 
large cirrus canopy in acivance of the frontal band 
creates a significant overesti nation of the precipitation 
pattern, projectin'’ the leading edge of precipiation 
substantially to the east. Orientation is also distorted 
since the IT enhancement slants northwest to southeast 
v;hile the radar depiction slants north to south. Cloud- 
top temperatures, ranging from 213 K to 243 2, imply 
changes in precipitation intensity throughout the region, 
hov/ever, all reported intensities are steady and light. 

As seen in Fig. 15, the h-PS technique is classified 
as fair agreement. The enhancement curves used in the 
model do not completely differentiate between the cirrus 
shield and the frontal band which causes the same over- 
estimation, but not as gross as stated in the IT 
technique. Area coverage is 25 T. Orientation of the area, 
unlike the previous technique, is more north to sout’n 
like the radar summary chart. Intensity over the observed 
area remains nearly constant v;hich is consistent v/ith the 
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steady licht precipitation pattern depicted on the radar 
summary chart and, therefore, is classified as correct. 
Yet the leading edge of the precipitation appears 
substantially to the east of the observed edge. 

The bispectral teclinique observed in Fig. 16 shoMS 
the sane problems indicated in the IR threshold 
technique. Again, area coverage is 20F, \/hich classifies 
tlie technique as having poor agreement. The cirrus shield 
causes the sane overestimation. Orientation of the 
bispectral precipitation pattern axis remains northwest 
to southeast \;hich differs from the observed pattern. 
Intensity is not distinguishable since there is little or 
no shading differences in the frontal zone or cirrus 
shield. The leading edge of the precipitation is again 
substantially east of the observed edge. 

Table 1 contains a summary of case number 1. All 
methods are marginally poor cases. The h'PS technique 
seems to do slightly better than the other three 
satellite estimates. I/here all of the other methods 
overestimate the amount of precipitation occurring, the 
h’ESDIS model severely underestimates it. Alignment is 
incorrect in all techniques except the MPS. Both the 
h'ESDIS and NPS methods show correct precipitation 
intensities while the IR technique overestimate its and 
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13 indicates that 



the bispectral underest inates it. Fi^. 
the XESDIS technique correctly identifies the leading 
edge of precipitation v/hile all of the other satellite 
estimates show the leading edge substantially to the 
east. 

C. SUMXAFY OF ALL CASES 

Table 13 summarizes the evaluations of all tv/elve 
cases (Tables 1-12) and provides an overall comparison of 
the techniques using the calculations derived from, all of 
the cases. Although the percentages do not show many 
dramatic differences betv;een the techniques, the XESDIS 
technique is shov/n to be superior in almost every 
category. The m.ost likely explanation for this improved 
capability is the human input in the analysis stategy. 

In percentage of area coverage, the IR threshold 
technique is the least effective (37%). But the other 
three techniques only shov; modest to good skill and have 
scores \/ithin ten percentile points of each other (i.e., 
52% to 61%). Table 13 shows that all three SPADS 
techniques consistently overforecast the amount of 
precipitation, some more than others. But the h'ESDIS 
technique, as seen in Table 13, over- and underestimates 
precipitation equally. Again, inserting human judgement 
into the analysis strategy seems to improve the overall 
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area estimation. Another explanation for the low scores 
seen in all of the techniques is that they nay be 
artifically lower than normal because the selected area 
of interest did not include the v/estern edge of 
precipitation. And it v;as over v/ater, along the western 
edge of the precipitation, v;here the technique seems to 
make its best analysis. Thus, given the area coverage 
scores alone, one would incorrectly conclude that none of 
the teciiniques provide better than good representation of 
the observed precipitation patterns. 

In addition. Table 13 shows that all techniques do 
an excellent job in getting the correct alignment, but 
only the NE3DIS and hPS methods show adequate skill in 
predicting intensities of precipitation. The results of 
the i'lPS model are extremely promising since there is no 
man involved in the analysis strategy. The IR technique 
normally overestimates v;hile the bispectral often 
underestimates the intensities. Intensities were 
extremely difficult to evaluate for the bispectral 
technique. The shades of yellov; and orange proved much 
more difficult to discern than shades of gray. Another 
factor was the quality of hardware. Poor resolution on 
the CRT or hardcopy print device caused further hardships 
in estimating precipitation intensity. 
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Finally, Table 13 indicates that tlie IJESDIS technique 
shows excellent skill in depictin^^ the leading edge of 
the observed precipitation. The SPADS satellite 
estimation techniques seen to consistently place the 
leading edge substantially to the east. 

Fig. 17 looks at the area coverage percentage 
scores for all techniques and cases over the v; inter 
season from Pec 19 8 5 - Tar 1937. Although no specific 
pattern is found v/ith this graph, it is easy to pick out 
the cases (e.g., 2, 4 and 5) v/here one particular 
technique did significantly better than the others. It 
also highlights a gradual improvenent of the SPADS 
techniques as the inter season progressed which is 
probably attributable to a strengthening of inter 
storms. Therefore the stronger storms were faster moving 
and more sharply defined than the earlier slow moving 
sy s t ems . 

Table 14 shows the TESDIS category used in each case, 
b'otice that only two (i.e.. Cloud Band and Subsynoptic- 
Scale h'ave) of the five basic NCSDIS categories mentioned 
in Chapter IV were used in all twelve cases. This 
probably stemmed from the initial criteria for selecting 
cases v/hich required precipitation on the radar summary 
chart somewhere within the radar umbrella of either 
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Sacramento, CA or Medford, 02 . To accomodate this 
requirement, several of the other categories, v/hich occur 
farther to the north, \;ere excluded. 

Fig. IS takes all the cloud-top temperature data 
collected from all tv/elve cases and plots it out in 
a graphical format. As observed in the figure, the top of 
the mean cloud-top temperature range is approximately 241 
v/hich is close, although somev/hat colder, to the 
threshold range (i.e., 246 to 253 IC) used in models, 
sighted by Paul (1983). The explanation for this 
difference is not obvious, but may be explained by the 
v/inter season and higher latitude v/here the cases were 
selected . 
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VI, co^XLrsio^:s atd n:-cn:::;7.:nATio:?s 

The follov;in 3 conclusions arc derived from results 
described in Chapter V: 1) the manual (man intensive) 

h^STIS technique proved superior to the objective SPADS 
techniques, 2) all techniques overestimated observed 
precipitation, but the ?;PSDIS method reduced th.e amount 
of overestimation and even underestimated by 
approximately the same amount, 3) all techniques did a 
better job representing the western edge rather than 
eastern edge of precipitation, A) all techniques shov; 
excellent alignment, 5) only the PESDIS and !IPS methods 
show sl:ill in predicting precipitation intensity, 6) the 
PESDIS technique shov/s excellent skill in depicting the 
leading edge of precipitaion, 7) the average range of 
cloud-top temperatures which match radar summary chart 
patterns is approximately 5^ C colder than previous 
studies indicated* 

Precipitation remains one of the most difficult 
\;eather parameters to monitor. Current conventional 
methods (surface and radar observations) of reporting 
precipitation provide excellent coverage over most of 
the United States. They provide a timely product \/hich 
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is used daily in business, transportation and recreation, 
II ov; ever, the coverage ends approximately 50 nautical 
miles offshore and do not aid the Ilavy’s needs to 
monitor precipitation events across the v/orld^s oceans, 
':"his thesis illustrates that using digital data from 
geostationary satellites, combined v;ith proven algorithms 
and man/ machine techniques, provides a reliable product 
v;!iich can be used in the fleet. 

The hhSDIS technique has been used operationally at 
several locations across the United States, Its 
performance is proven and regarded as a technique v;orthy 
of transmission over the UUIFAX circuit. Therefore, it is 
r ecom*m;ended that a thorough study be made of this 
technique, evaluating all five categories. The evaluation 
should be made over the open ocean in addition to an 
overland controlled region. If the evaluation proves that 
the technique provides excellent correlation with 
observed precipitation patterns, it is further 
recom. mended that a technical manual be developed for use 
by the '’aval Oceanography community. The manual would be 
available to assist Oceanography centers in producing a 
precipitation analysis for transmission over the fleet 
broadcast . 
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TABU': 1 

lonipa r i son of Tecliniqiios for Case Bumh 
at 1830Z on 4 December 1987 
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TABLE 3 

Comparison of Toclinifiuos for (!ase [lumber 
at 1830Z on 2(3 December 1937 
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TARLK 4 

Couparison of Tecliniques for Case Hunber 
at 1830Z on 1 January 1987 
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TAl’.LI', 5 

Compariaon of Techniciiies for f!nao Ilumbor 
at 1030Z on 2 January 1987 
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TABLE 6 

Comparison of Tecliniques for Case Iluinl)er 6 
at 1830Z on 3 January 1087 
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TABLE 7 

Comparison of Techniques for Case Humber 
at 1B30Z on 22 January 1987 
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TABLE 8 

Comparison of Tecliniques for Case Number 
at 1830 Z on 2 February 1987 
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TABLE 9 

Conp arisen of Tecliniqiios for Case Number 
at 213 OZ on 9 February 1987 
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TABLE 10 

Comparison of Tecliniques for Case Number 10 
at 213 OZ on 12 February 1907 
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TABLE 11 

Comparison of Techniques for Case !Ium])er 11 
at 1830Z on 13 I'ebruary 1087 
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Comparison of Techniques for Case i!umber 12 
at 0030Z on 13 March 1987 
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TABLI' 13 

Summary of Tables 
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TABLE 14 

NESDIS Cater^ory and Subjective Skill Compared 
with Radar and Surface Observations 
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Figure 1. GOES Visible Mercator Image for 
2130Z on 9 February 1987 
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Figure 2. GOES Infrared Mercator Image for 
2130Z on 9 February 1987 
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Figure 3. ’3IC Radar Summary Chart for 2135Z 
on 9 February 1907 
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Figure 4. 



Infrared Schematic of the MESDIS Cloud Band 
Catcf’ory from Scofield and Spayd (1934) 
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SUB-SYNOPTIC SCALE WAVE CATEGORY 

IN «N OVERRUNNING ZONE 
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Figure 5. Infrared Schematic of HESDIS Subsynoptic- 
Scale '.’ave Category in Overrunning Zone 
from Scofield and Spayd (1984) 
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Figure 6. ^'ESDIS Technique overlayed on High 
Resolution GOES-6 Visible Image for 
2100Z on 9 February 1987 
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Figure 7. Infrared Technique (-25*^ C to -55*^ C) 

Image for 213QZ on 9 February 1987 
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Figure 8. NTS Model Results for 2130Z 
on 9 February 1987 
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Figure 9. Bispectral Image for 2130Z 
on 9 February 1987 
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Figure 10. GOES Visible Mercator Image for 
1S30Z on 4 December 1986 
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Figure 11. GOES Infrared Mercator Image for 
1S30Z on 4 December 1986 
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Figure 12. IFIC Radar Suminary Chart for 1S30Z 
on 4 Decenber 1980 
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Figure 13. NESDIS Technique overlayed on High 
Resolution GOES-6 Visible Inage for 
1901Z on 4 December 19G6 
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Figure 14. Infrared Technique (-25° C to -55° C) 
Image for 18302 on 4 December 1986 
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Figure 15. NFS Model Results for 1S30Z 
on 4 December 1987 
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Figure 16. Bispectral Image for 1830Z 
on 4 December 1986 
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Fioiire 17. Comparison of Area Coverage by 
Techniq ue 
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Figure 18 . Cloud-fop Temperatures Correlating 
v/ith Observed Precipitation 
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